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Abstract 

A general condition for sharp transition of decay rate from quantum to 
thermal regimes is derived in dissipative tunneling models when position- 
dependent mass is involved. It is shown that the effect of dissipation in general 
changes the order of the phase transition. Especially, for the models with con- 
stant mass the Ohmic dissipation enlarges the range of parameters for first or- 
der phase transitions. In the case of second-order phase transition the Ohmic 
dissipation suppresses the decay rate near the transition temperature(Tc). For 
the super-Ohmic case the dissipation yields an opposite effects to the Ohmic 
dissipation within exponential approximation. 
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I. INTRODUCTION 



It is well known that the decay of the metastable states at zero temperature is deter- 
mined by pure quantum tunneling process whose dynamics is described by classical configu- 
rations called bounce in Euclidean space . Since, however, the decay rate is determined at 
high temperature by thermal activation which corresponds to classical configuration called 
sphaleron [Q, there exists some temperature at which the transition from classical- to 
quantum-dominated decay occurs. This phase transition problem was firstly discussed by 
Affleck within quantum mechanics. He demonstrated that under certain assumptions 
for the shape of the barrier the transition between the thermal and quantum regimes is 
dominated by solutions with a finite period in the Euclidean time that smoothly interpo- 
late between the zero-temperature bounce and the static high-temperature sphaleron. The 
transition is thus a second-order one. Chudnovsky [|], however, has shown that the type of 
the phase transition in the crossover from thermal activation to thermally assisted quantum 
tunneling is completely dependent on the shape of the potential barrier. He also has shown 
that the order of the phase transition is easily conjectured by P-vs-E graph, where P and E 
are Euclidean period and energy, repectively. The sharp crossover between the thermal and 
thermally assisted tunneling occurs when P = P{E)-cnTve possesses a minimum at E = Ec, 
which is different from the energy of sphaleron solution. Based on Chudnovsky's observation 
the sharp first-order transitions are found at spin tunneling systems with and without 
external magnetic field. 

Recently, a sufficient criterion for the first-order phase transition in decay problem of the 
metastable state is obtained by carrying out the nonlinear perturbation near the sphaleron 
solution in the two-dimensional string model 0. Inspired by spin-tunneling problem, the 
result of Ref. 0] is subsquently extended to the quantum mechanical model when mass is 
position-dependent 0. 

The purpose of this paper is to derive a general criterion for the sharp first-order phase 
transition in the quantum mechanical tunneling models when the position-dependent mass 
and dissipation are involved. We first consider a system with Ohmic dissipation in which the 
spectral density J{uj) is linearly proportional to frequency, i.e., J{uj) = au where a is the 
dissipation coefficient. In this case Caldeira and Leggett already derived the effective action 
in their seminal paper |T0] . We will use their result for the derivation of the criterion. The 



extension to the case of super-Ohmic dissipation where J{oj) = au^ is also examined. In 
the super-Ohmic case a counter term corresponding to the deformation of potential due to 
environment is introduced. The absence of the counter term in the Ohmic case is due to our 
use of the result of Ref. 0], in which a proper counter term is already used for the Ohmic 
dissipation. Since we work at Euclidean space, we use a prescription for the dissipation 
term which does not break the time-reversal symmetry in the both Ohmic and super-Ohmic 
cases. Carrying out the nonlinear perturbation we will derive the general condition for 
the first-order phase transition in the dissipative quantum mechanical system when mass is 
position-dependent in Sec. II. It is found that the effect of dissipation is only deformation of 
the eigenvalues for the temporal fiuctuation operator / defined in this section. To get some 
physical intuition we will apply this general criterion to the simple quantum mechanical 
models in Sec. HI. Brief conclusion and the direction of the future research in this field will 
be given at the final section. 
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II. CRITERION OF FIRST-ORDER PHASE TRANSITION FOR DISSIPATIVE 

QUANTUM SYSTEMS 



The effect of dissipation on quantum tunneling is investigated in Ref. [|T0| by introducing 
the infinite number of oscillators as an environment and assuming the linear coupling between 
the environment and system. We follow their formalism to obtain the effective action for 
the dissipative system. 

The effective action in Euclidean space is given as 

S[q{r)] = / [-M{q)f + V{q) + 6V{q)]dT + - / dr' dr^{r - r'){q{r) - q{r')}\ (1) 
Jo 2 / J-oQ Jo 

where 

^{r-r')^^ r J{u:)e-\--^'\du. (2) 

Here, the mass of quantum system is generally taken as position-dependent which is moti- 
vated from spin tunneling models The final non-local term in Eq.(|l|) represents the 
effect of dissipation. It is emphasizing to note that the counter term 5V{q) is introduced in 
Eq.(P |ll|]. In fact, the above effective action without 5V{q) has been derived in considera- 
tion of a counter term which cancels divergence in the case of Ohmic dissipation. However, 
for the cases of non-Ohmic dissipation new divergence is expected to appear, so that appro- 
priate confer term may have to be needed. Shortly, it will be shown that this is a generic 
case. 

Now, we take the spectral density generally as 

J{uj) = au"" (3) 
where n is an positive integer. Then, the effective action can be written as 

where r(n) is the Gamma function. From this effective action, the equation of motion is 
given by 

Z TT J-oo \ T — T' I"^"-^ 

where the prime denotes the derivative with respect to q. In the following subsections we 
will derive the criterion of the first-order phase transition in the cases of Ohmic dissipation 
and super-Ohmic dissipation ( J (u;) = auj^). 



A. The Ohmic dissipation 

In the case of Ohmic dissipation, the spectral density is linearly proportional to frequency. 
The equation of motion is, thus, Eq.(S) with n = 1. The dissipation term can be rewritten 
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in terms of the Fourier transform partner of 5'(t), i.e., q{uj), which makes the equation of 
motion to be 

M{q)q + l^^q' - « duqiu) \ u \ e'^^ = V\q). (6) 

In the derivation of this equation we take a prescription for the dissipation term 

_ a ^^, q{T) - q{T') _^ _a qjr) - qjr') 

71 J~oo (r - r')2 TT i-oo (r - r' + ie)(r - r' - ie) ' ^ ^ 

which preserves the time-reversal symmetry, and 6V' = as mentioned before. Since, in fact, 
the dissipation breaks the time-reversal symmetry in real space-time, another prescription 
was used at Ref. [|lT|. However, in Eclidean space the dissipation does not yield damping 



motion, which justifies our choice of the prescription. 

Now, following Ref. , let us determine the type of transition by expanding the equation 
of motion at the sphaleron solution as 

(li'r) = qs + aT]i{T) (8) 

or equivalently 

q{u) = qs + afii{uj), (9) 

where a represents an oscillation amplitude near sphaleron solution. Since it is sufficient in 
determining the order of phase transition to consider only the solutions near sphaleron, we 
can assume a is very small constant. Substuting these expressions into Eq.(^ the equation 
of motion within the first order of a becomes 

(/-%i(r) = 0, (10) 

where operators / and h are defined as 

/ = M(g,)^ -^j du\u\ j dTe-'^\ 

h = V"{qs), (11) 
respectively. In order to solve this equation, we take a trial solution 

ryi(r) = cosfir (12) 

and equivalently 

fl^{uj) = ]^{5{Sl-uj) + 5{Sl + uj)). (13) 

Using this trial solution the frequency VL near the sphaleron solution is obtained within the 
first order of the amplitude a as 



where 

= -MM- 

Comparing | Qq^ \ with Ug which is a ^ hmit of I n'o^ I, it is easy to understand that the 
Ohmic dissipation reduces the frequency. Since in the case of second-order phase transition 
the transition temperature is determined by this | Qq^ |, i.e., =| Qq^ \ /27t, this result 
imphes that the Omic dissipation decreases the transition temperature. As shown in Fig.l, 
the decrease of Tc means increase of action value and, then, suppression of decay rate near 
Tc. We will, however, show that the super-Ohmic dissipation gives the opposite effect on 
decay rate within exponential approximation. 

Next order calculation can be conducted by taking 

q{r) = qs + ar]i{T) + a'^r]2ir). (16) 

Then, we find equation of motion for ?72(t) as 

a{i-h)7]2{T) = -{l-h)r],{T) + aW,{T), (17) 

where 

Wi{t) = {^fM'{qs) + ]^V"'{qs)) cos^ fir - ^fi^M'(g,) sin^ fir. (18) 

Since (/ — h) is an hermitian operator, taking a scalar product with | ?7i > on both sides of 
the equation of motion yields 

a(/(fi) - h{Q)) <% \r]2 >= - < % \ 1 - h \ t]i > +a < % \ Wi >, (19) 

where 

/(fi) = -n^M{q,)-a I fi I, 

h{Q) = V'iqs). (20) 

As usual perturbation theory, < r/i | 772 > is zero. Furthermore, the second term of the right- 
hand side in Eq. ([T9|) is zero because of r-integration. Then, the above equation(Eq.([T9|)) is 
identical with Eq.(^). Therefore, within the second order of a we can not find the variation 
of the frequency from fio\ i.e., 

fig) = fig). (21) 

Before calculating next order frequency fi^ , we would like to evaluate ?72(t) exphcitly for 
later use. This is achieved directly from Eq.(|T^, which yields 



V2{r) = ii-h)-'WiiT) 

= 91+92 cos 2fir, (22) 



where 
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n^M'{q,) + V"'{q,) 

^ 3n'M'{q.,) + V"'{q,) 

A[AM{qs)Q^ + a\2Q\+V"{qs)y ^ ^ 

Now, taking a third order correction into g(r) as 

= qs + ar]i{T) + a^r]2{T) + a%{T), (24) 

and inserting the above expression into Eq.(^), the equation of motion for 773(1") is straight- 
forwardly obtained as 

a\l-h)r]3 = -{i'h)r], - a{l - h)r]2 + aW^ir) + a^W2{r), (25) 

where 

W2{t) = {Q^g,M'{qs)+giV"'{q,))cosQT 

+ {5n^g2M'{qs) + g2V"' (qs)) cos 2nT cos 

+ hn^M"{q,) + ]-V"") cos' fir (26) 
-2Q'^g2M' (gj sin Qt sin 2Qt 
~Q'^M"{q,) cos Qt sin^ Qt. 

Scalar product with r/i in Eq.(^) yields an equation 

- {1{Q) - h{Q)) <r]i\Vi> < r^i | >= (27) 

(3) 

in which we can determine Qq . Since the first-order phase transition between thermal and 
thermally assisted quantum tunneling regimes occurs when the period of solution near the 
sphaleron increases with approaching to the sphaleron solution, we get a condition for sharp 
transition by | Qq} \ > \ Qq \ . This is identical to the condition that the value of the left- 
hand side of Eq.(|3) at fi =| Q^q^ \ is negative. Therefore, we finally obtain the criterion of 
first-order phase transition as 

n'i^'ligi + l92)M\q,) + ^M"(g,)] + {g^ + ^^?2)V^'"(gs) + ^V""(g.) < 0, (28) 

where gi and g2 are evaluated at Q = Qq\ Since the dissipative coefficient is involved at 
Q^^\ gi, and g2, the criterion of first-order phase transition is generally different from non- 
dissipative case. Particularly, when particle mass is constant, the dissipation coefficient a 
appears only in g2, and from the fact V"{qs) < we can conclude that the Ohmic dissipation 
enhances the possibility for the occurrence of the sharp first-order transition. 



B. The super-Ohmic dissipation 

In this subsection we consider the criterion for the first-order phase transition when the 
dissipation is super-Ohmic, i.e., J{uj) = au', which appears in the analysis of macroscopic 
magnetization tunneling fl^ . 
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The dissipation term in Eq.(|^) with n = 3 can be written again in terms of q{uj). After 
integrating it at complex r'-plane by using same prescription, one can find the equation of 
motion (Eq.(^)) to be in the form 

Miq)q + l^^^q^ + a T dujq{u) \ u 1=^ e'^^ = V\q). (29) 
Z oq J —oo 

In deriving Eq.(|29|) we used 5V' = ^q to cancel a divergence which appears in the course 
of integration. The presence and absence of divergence in super-Ohmic and Ohmic cases, 
respectively, are because that we start with the effective action (Eq.(|l|)) as mentioned in the 
previous section. 

The remaining calculation for super-Ohmic case is equivalent to that for the Ohmic case. 
The differences are foUowings. The / operator is changed in this case into 

/ = M{qs)^ + ^lduj\u\' I dTe-'^\ (30) 

and the frequency near the sphaleron solution within the first order of the amplitude a, fi^ ^ , 
becomes the root of the equation 

M{q)^f -a\n\^ +V" = 0. (31) 

This equation tells that for the second-order phase transition the transition temperature 
Tc =1 ^l^g^ I /27r becomes higher than that of the non-dissipative case, uJs/27t. This means 
that contrary to the Ohmic dissipation case, the tunneling rate near is enhanced by the 
super-Ohmic dissipation. 

Next order calculation shows that r]2 has a same form with the previous case(Eq.(P^), 
where gi is given at Eq.(^) and g2 is changed in the form 

^ Sn'M'jqs) + V"'{qs) , . 

4[4M(g,)^]2 _ a I 2fi |3 +V"{qs)] ' ^ ' 

Finally, the criterion of first-order phase transition is equivalent to Eq.(ESD except that g2 



is replaced by Eq.(P^ and Eq.(PB|) should be evaluated aX Vt = q!^^ instead of ■ 



then, evident that when particle mass is constant the super-Ohmic dissipation reduces the 
possibility for the occurrence of the sharp first-order phase transition within exponential 
approximation. 



III. APPLICATION TO QUANTUM MECHANICAL TUNNELING MODELS 

A. Asymmetric double well case 

Consider a usual double well potential with a symmetry-breaking term proportional to 
i.e., 

V{q) = \{q'-lf-fq\ (33) 
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This type of potential is frequently used in field theories for the application to cosmology 
In this case the sphaleron solution is 



qs = 0. (34) 

We assume that particle mass M is constant, so that the derivaltives of mass are zero. 
Substuting the derivatives of the potential at qs 

V"{qs) = -2, V"\qs) = -6f, V^""(g.) = 12 (35) 

and Eq.(^) into Eq . (|28|) , it is easy to obtain following criterion 



2^Vm' + M ~ > 2 + 2(3/2 + 1) ^^^^ 

for the first-order phase transition when the dissipation is Ohmic. Unfortunately, for any val- 
ues of positive a and M the left-hand side cannot exceed 2, i.e., always second order. Hence, 
dissipation does not change the order of the phase transition. This result is maintained even 
in the super-Ohmic case. 

We can also obtain the criterion for a double well potential with q asymmetric term as 

V{q) = ^{q'-ir-Fq. (37) 

In this case, the result is very similar to q^ case. There is no first-order phase transition. 
This result can be expected since the both asymmetric potentials Eqs.(^) and (^7\), when 
the dissipation term is ignored, give actions proportional to each other, which can be easily 
shown by a appropriate translation of q and scahngs of q and r. 



B. Case for spin tunneling-inspired model 

Consider an hamiltonian 

where 

M(0) = \ ■ , .. , (39) 

2Ai(1 — Asm (p) 

and 

V{(l)) = K2S{S + 1) sin^ (j). (40) 

Although this hamiltonian can be derived by the coherent representation from the hamilto- 
nian of spin tunneling model 

H = K^Sl + K2SI (41) 



where Ki and K2 represent an anisotropic constants and \ = K2/ Ki, we will not go through 
the content of physics on spin tunneling in this paper. Instead, we will restrict ourselves 
into the discussion on the effect of dissipation in the hamiltonian, Eg. (|38D . 
Since the sphaleron solution is simply 



TT 

2' 



(42) 



it is easy to obtain 



V"{<f)s) = -2K2S{S + 1), V"'{<j)s) = 0, V""{<j)s) = 8K2S{S + 1), 

MUs) = 7^ r, M'(d),) = 0, M"(0J = — —- — (43) 

Hence the use of Eq. (pSD for the Ohmic dissipation yields the condition for first-order phase 
transition as 

S{S + 1) 1-^ 



This result is shown in Fig. 2. It is shown that, unlike the case of constant mass, the effect of 
the Ohmic dissipation reduces the range of parameters for the first-order phase transition. 

Since it is impossible to get an analytic expression of criterion for the super-Ohmic case, 
one has to resort to the numerical calculation. The type of transition in the parameter space 
is given at Fig.3. The effect of the super-Ohmic dissipation enlarges the parameter range 
for the first-order phase transition, which is opposite behavior to the Ohmic case. 

It is worthwhile noting that the above criterion may not be realized in real spin system. In 
order to study the effect of dissipation in the real spin system, one has to derive the effective 
hamiltonian from the original hamiltonian(Eq. (^ID ) with the dissipation term fully through 
coherent representation or particle mapping |jl5|. Since this procedure in general produces 
a complicated dissipation term, our cirterions does not hold in this case. Therefore, what 
we have done in this subsection is to obtain the general criterion for the sharp transition for 
the model hamiltonian (Eq.(^8|)) with Ohmic and super-Ohmic dissipation and to show the 
possibihty of the change in type of phase tansition due to dissipation. 



IV. CONCLUSION 

Performing the nonlinear perturbation, the general condition for the sharp first-order 
phase transition of decay rate between thermal regime and thermally assisted quantum tun- 
neling regime has been derived when position-dependent mass and dissipation are involved. 
It has been shown that in the models with constant mass Ohmic dissipation enhances the 
possibility for the occurrence of the sharp transition, while super-Ohmic dissipation reduces 
it. Application to two simple quantum mechanical models is given at the previous section. 
In the asymmetric double well case the phase transition is always second order regardless 
with or without dissipation. In this case, by comparing ^Iq^ and 0,^'^ with a;^ one can per- 
ceive that the Ohmic and super-Ohmic dissipation suppresses and enhances the decay rate, 
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respectively, near the transition temperature Tc within exponential approximation. Similar 



feature at zero temperature has been reported in Ref. ||TT . 

For the case of J{oo) = aio^ it is not clear how to give a proper prescription which does 
not break the time-reversal symmetry. This difficulty does not seem to originate from the 
property of Euclidean space. In real time space, a retarded Green function is generally 
choosed in order to break the time-reversal symmetry. However, when J{lS) = au"^, in spite 
of using the retarded Green function, the time-reversal symmetry is not broken. 

The generalization of our result. Eg. (pHf) , to nonlinear coupling between an environment 
and system and general type of dissipation might be highly non-trivial. In this case the 
dissipation term in the equation of motion, Eq.(|^), is proportional to 



TT J-oo \ T — t' 



where Zi and Z2 are some constants. Hence, the integration with respect to r' through Fourier 
transform as was done in Eq.(^) is impossible. In our opinion the only way to break through 
this difficulty is to rely on the numerical analysis. This will be discussed elsewhere. It is also 
interesting to apply the present method to the spin tunneling problems. Since dissipation 
arises from the coupling with an environment which contains quasi-particles like photons 
or phonons, it is possible to construct a hamiltonian which contains the spin-environment 
coupling. Then, by integrating over the environmental degrees of freedom and with the help 
of coherent representation or particle mapping, one can derive an effective hamiltonian in 
which the effect of dissipation is fully involved as a non-local term. From this procedure, 
one can investigate the effect of the dissipation on the phase transition in spin tunneling 
problem. This will also be discussed elsewhere. 
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FIGURES 



FIG. 1. The seond-order phase transition. The decrease of transition temperature means the 
increase of action value near Tc, indicating suppression of the decay rate. 

FIG. 2. The parameter domains for first and second-order transitions in the hamiltonian 
(Eq.(40)) with Ohmic dissipation. 

FIG. 3. The parameter domains for first and second-order transitions in the hamiltonian 
(Eq.(40)) with super-Ohmic dissipation. 
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